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ABSTRACT 

The hydrogen-deficient, carbon-rich R Coronae Borealis (RCB) stars are known for being prolific 
producers of dust which causes their large iconic declines in brightness. Several RCB stars, including 
R CrB, itself, have large extended dust shells seen in the far-infrared. The origin of these shells is 
uncertain but they may give us clues to the evolution of the RCB stars. The shells could form in 
three possible ways. 1) they are fossil Planetary Nebula (PN) shells, which would exist if RCB stars 
are the result of a final, helium-shell flash, 2) they are material left over from a white-dwarf merger 
event which formed the RCB stars, or 3) they are material lost from the star during the RCB phase. 
Arecibo 21-crn observations establish an upper limit on the column density of H I in the R CrB shell 
implying a maximum shell mass of <0.3 M 0 . A low-mass fossil PN shell is still a possible source of 
the shell although it may not contain enough dust. The mass of gas lost during a white-dwarf merger 
event will not condense enough dust to produce the observed shell, assuming a reasonable gas-to-dust 
ratio. The third scenario where the shell around R CrB has been produced during the star’s RCB 
phase seems most likely to produce the observed mass of dust and the observed size of the shell. But 
this means that R CrB has been in its RCB phase for ~10 4 yr. 


1. INTRODUCTION 


R Coronae Borealis (R CrB) is the prototype for its 
eponymous class of stars, which are very rare and have 
many unusual characteristics including extreme hydro- 


of 8 magnitudes or more (Lambert & Rao 

1994 

Clay- 

ton 

1996 

2012 

). These declines are caused 

by clouds of 


which are later dissipated by radiation pressure. Only 
about 100 R Coronae Borealis (RCB) stars are known in 
the Galaxy. Therefore, these stars may be the result of 
a rare form of stellar evolution or are in an evolutionary 
phase that lasts only a short time. 

Two scenarios have been suggested for producing an 
RCB star: the double degenera te (DD) and the final 


helium-shell fla sh (FF) models (Iben et al. 1996 Saio 


& Jeffery 2002). In the DD model, an RCB star is the 


result of merger betw een a CO- and a He-white dwarf 
(WD) ( |Webbink||1984[ ). The derivation of CNO isotopic 


vors the DD scenario for these stars i 

Clayton et al.|2005 

2007 Garcfa-Hernandez et al. 

2009 

2010|). In the FI 


final helium-she ll flash and expands to supergiant size 
In this scenario, the star may have 
irough a planetary nebula (PN) phase. 


(Fujimoto 1977). 
recently gone t 


Three stars (Sakurai’s Object, V605 Aql, and FG Sge) 
have been observed to undergo FF outbursts that trans¬ 
formed them from hot evolved stars into cool gi ants w ith 
spectroscopic properties similar to RCB stars (Clayton 
fe De Marco |1997l_ Gonzalez et al. |1998| Asplund et al. 
1999[ 20001 (Clayton et al. 2006). These FF stars 


1998 


are al 


surrounded by Pi\e. 
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Several RCB stars have extended dust shells that are 
seen in reflecte d light in the visible or in emission from 
the far-IR (e.g.,|Schaefer|1986||Walker|1985[|1986||Bright| 
et al.|2011|). Infrared spectroscopy of RGB stars has been 
bl 


possible thanks to the Infrared Space Observatory (ISO) 
and the Spitzer Space Telescope, which has permitted 
the extraction of the characteristics of the IR emitting 
dust shell around RCBs (Lambert et al. 2001; Garcfa- 
Hernandez, Rao & Lambert 2011a,b, 2013). It has been 
suggeste d that the mass loss in R CrB and V854 Cen is 
bipolar ( Rao &; Lambert|1993 Clayton e t al. 1 199 71) . Thi s 
geometry is very common in PJNIe (Baiick & Frank 2002). 
Further, the nebulosity, including cometary knots, seen 
around R C rB and UW Cen has similar morphology to 

Clayton et al. 


2011 


2015, 


a PN shell (Clayton et al. 
in preparation). In this paper, we investigate whether 
the shell around R CrB is consistent with the scenarios 
suggested for the evolution of the RCB stars. 

2. ARECIBO 21-CM OBSERVATIONS 

Radio observations containing the region around 
R CrB have been obtained from the first data re¬ 
lease (DRlrl) Galacti c Arecibo L-band Feed Array H I 
(GALFA-HI) Survey ([Pee k et al. fMl] ). The GALFA-HI 
survey provides both high resolution (4 7 ) and high sensi¬ 
tivity (typical rms ~ 80 mK) due to their 305-m aperture 
and the installation of th e Arecibo L-band Feed Array 
(ALFA) ( Peek et ah|2011 ). Fully reduced “narrow” and 
“wide” band data cubes were retrieved from the survey. 
The former provided 0.18 km s _1 resolution in the local 
standard of rest (LSR) in the range vlsr = ±190 km 
s _1 , while the latter provided 0.74 km s _1 resolution in 
the range vlsr = ±750 km s _1 . A full description of 
the data acquisition, which includes both drift and bas- 
ketweave scanning, and reductions f or th e GALFA-HI 
survey are presented by (Peek et al.| |2011) (see their §3 
and §4, respectively). 

3 https://purcell.ssl.berkeley.edu/ 























































































2 


3. ORIGIN OF THE R CRB DUST SHELL 

3.1. The R CrB Dust Shell 

The foreground extinction toward R CrB is quite small, 
E(B-V)^0.035 m ag, since it lies at high Gal actic lat itude 
(b J/ =±51°) (Schlegel et ah 1998 Schlafly & Finkbeiner 
20111. At ma ximum light, R CrB is V=5.8 mag and 
B-V =0.6 mag (Lawson et al. 1990). Based on the abso¬ 
lute magnitude/ettective temperature relationship found 
for the Large Magellanic Cloud RCB stars, the abso- 
lute magnitude of R CrB is estimated to b e My=-5 mag 
(Alcock et al.|2001 Tisserand et al.|[2009). For the anal¬ 
ysis in this paper, we thus adopt a distance to R CrB 
of 1.4 kpc. The large extended far-IR shell, with a ra¬ 
dius of ~10', surrounding R CrB was discovered with 
IRAS and then studied further with the Spitzer and Her- 
schel telescopes (Gillett et al. 1986; Clayton et al. 2011; 
Garcfa-Hernandez, Rao & Lambert 2011b, 2013). At the 
assumed distance of R CrB, the radius of the shell corre¬ 
sponds to 4 pc. Monte Carlo radiative transfer modeling 
of the R CrB shell s ugges ts that it contains 10~ 2 M 0 of 
dust ( Clayton et al.pOlT ). 


3.2. A Planetary Nebula Shell? 

The PNe around the FF stars, V605 Aql, Sakurai’s 
Object, and FG Sge are still ionized. The shell around 
R CrB is not. There is a small subclass of RCB stars 
that are much hotter (T e //= 15,000-20,000 K) than the 
typical RCB stars and are surround ed by PNe (Pollacco 


et al. 


1991 De Marco et al. 2002). They are not hot 


enough at present to ionize their surrounding PNe, but 
the nebulae have not had time to recombine. The shells 
around the cooler RCB stars could be old PNe that have 
recombined. There are some strong similarities between 
the morphology of the shells of UW Cen and R CrB, and 
some PNe such as the Eskimo Nebula. Cometary fea¬ 
tures seen in PNe such as the Eski mo are s imi lar to those 
seen in these two RCB stars (Clayton et al. :2011). The 
recombination time depends on the electron densities of 
the shells and could range from hundreds to thousands 
of years. If n e =200 cm" , as assumed for Sakurai’s Ob¬ 
ject, the recom binati on time for R CrB’s PN shell would 
be ~380 yr (Pollacco 1999). The amount of time that 
R CrB has been an RCB star is unknow n but it is one of 
the first va riable stars to be discovered (|Pigott & Engle- 
field 1797|), so we know that it has been in its RCB star 
phase tor at least 200 yr. 

The PNe seen around the FF stars, Sakurai’s Object 
and V605 Aql, are thought to be ~2x 10 4 yr old, both 
with expansion velocities of ~30 km s _1 and radii of 0.35 
pc for V605 Aql and 0.7 pc for Sakurai’s Object dPollacco 


et al. 1992 Guerrero & Manchado||1996 Pollacco 1999 b 
A FF star may not be old enough to produce a PIN shell 
of the size seen around R CrB as it only takes a few 
years f or the star to reach its RCB-like phase after the 
FF (pHerwig 2001). 

If the R CrB sh ell is, in fact, a PN shell then it should 
be hydrogen rich (fGillett et al. 1986). No pointed 21-cm 
measurements of any RCB circumstellar shells have been 
published previously. A visual inspection of the narrow 
band GALFA-HI data cube through velocity space was 
performed to see if there were any obvious H I features. 
The search focused on ±30 km s _1 of the vlsr of R CrB 
~ 37 km s _1 (The barycentric radial velocity of R CrB 



Fig. 1. - GALFA-HI field (2?25 x 1? 15) showing the location of 
the R CrB shell (green square) and the area used for background 
subtraction (blue square). The squares are 20' x 20 '. 


22.3 km s~ 4 (Lawson & Cottrell 1997)). No obvious 
structure was discerned. 

Next, a 20' x 20' square region, roughly the size of the 
dust shell, centered on R CrB was selected, as well as a 
background region of the same size as shown in Figure 
1. The background reference spectrum were subtracted 
from the spectrum toward the R CrB shell to determine 
if there is any H I emission that might be associated 
with the shell. While fluctuations in the background H 
I limit the comparison in many directions, one reference 
position, shown in Figure 1, has an H I spectrum that 
is identical within the noise to the spectrum averaged 
over the shell at the relevant velocities, and the difference 
spectrum (Figure 2) can be used to limit the H I mass of 
R CrB. 

The rms noise (cr rms ) in the 0.18 km s _1 channels of 
Figure 2 is 27 mK. Over a ±30 km s _1 interval, this gives 
a 3tr limit on column density, N# (atoms cm -2 ) in the 
R CrB shell of 4.9 x 10 17 . This was derived using N h = 
3.0 x 1.82 x 10 18 x u rms x dv x ^/n c h, where dv is the 21cm 
channel spacing (0.18 km s~ 4 ) and n c h is the nu mber of 
channels (^330) o ver the relevant velocity range (Dickey 
fe Lockman|1990 ). The mass of any neutral hydrogen, m 
solar masses, is then less than M# = 6.9 x 10" 28 x N h 
x D 2 x D, where D is the distance in pc (1400 pc), and 
D is the solid angle of the region in arcmin 2 . With our 
upper limit on N h, the corresponding 3er limit on H I in 
the R CrB shell is 0.3 M 0 . Assuming that the mass loss 
is contained in a shell that has been moving outward at 
30 km s _1 , the PN shell would take ~10 5 yr to expand 
to r= 4 pc. A shell of that age would almost certainly 
have had time to recombine and become neutral. 

The number density of H in this shell would be low, 
similar to that seen in the diffuse interstellar medium. 
Therefore, the fraction of H 2 will be near zero, and all o f 
the H will be atomic and neutral ( Rachford et al.|2 009). 
The 21-cm measurements presented here put a strong 
upper limit on the H mass in the R CrB shell. A study 
of the gas and dust masses in the shell of the PN, NGC 
6781, gives a total shell mass of 0.86 M 0 and a dust 
mass of 4 x 1 0~ 3 M^ and th erefore a gas-to-dust ratio 
of about 215 (Ueta et al. 20141. Assuming a gas-to-dust 
ratio of 200, the amount of dust in the R CrB shell would 
be <KT 3 M 0 . 


3.3. Mass Loss from a White Dwarf Merger? 

We have simulated the RCB DD formation scenario 
by running three binary WD merger models using our 
fully three-dimensional adaptive mesh refinement (AMR) 
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TABLE 1 
Model Results 


Parameter 

Model 1 

Model 2 

Model 3 

Mass Ratio (q) 

0.51 

0.70 

0.80 

Accretor Mass (M©) 

0.503 

0.526 

0.561 

Donor Mass (M©) 

0.255 

0.366 

0.45 

Total Mass (M©) 

0.758 

0.892 

1.011 

% Mass Lost 

2.3 

0.90 

3.3 

Initial Period (s) 

181 

118 

90.9 

Final/Initial Period 

9.9 

2.9 

6.9 

Initial separation (R©) 

6.30 x lO" 2 

5.00 x lO" 2 

4.38 x 10 -2 

Grid Box Dim. (R©) 

1.74 

1 

1.38 

Orbits driven at 1% 

2.5 

4 

2.5 

dx (R 0 ) 

5.66 x 10" 4 

6.51 x 10“ 4 

4.49 x 10“ 4 

Init. Sep. in Cells 

111 

77 

98 


Fig. 2.— A background-subtracted temperature vs. radial ve¬ 
locity plot for area around R CrB. The vertical green dashed line 
is the position of R CrB in radial velocity and the two vertical red 
lines mark ±20 and ±30 km s —1 from the velocity of R CrB. The 
feature at -70 km s —1 is part of a large elongated emission not 
associated with the R CrB shell. The horizontal green dashed line 
is at temperature = 0°. 


code (Marcello & Tohline, in prep). This code evolves 
density, total energy, and angular, vertical, and radial 
momenta on a rotating Cartesian mesh. Once started, 
the model is driven into closer contact by artificially re¬ 
moving angular momentum from the system at a rate 
of 1% per orbit for several orbits. Through AMR, we 
are able to run the models in large grids whose box di¬ 
mensions are 20-30 times larger than the orbital separa¬ 
tion. This allows for a more accurate determination of 
mass loss from the system than on a smaller grid, where 
much of the material exiting the grid does not possess 
enough energy to escape. The initial conditions were 
created using the self-consistent field technique (SCF) to 
create a detached synchronously rotating binar y WD in 
equilibrium just shy of beco ming semi-detached (Hachisu 
1986| |Even & Tohline|2009|. Each model run follows the 
merger of a (JO- and a He-WD in a close binary orbit. 
The q=0.8 and q=0.5 models have a h ybrid CO/He ac- 
cretor, which as described by Staff et al. (2012) is a WD 
with a CO core and a He envelope of about 0.1 M 0 . The 
results, along with the initial conditions, of our simula¬ 
tions are summarized in Table 1. 

In the models presented here between 0.9 and 3.3% of 
the initial mass escapes and is not expected to fall back 
onto the new merged star. Therefore, in the three models 
described in Table 1 the mass lost from the two WDs, 
which could form into a circumstellar shell around the 
new RCB star, is ~10 -2 M 0 . The gas that has escape 
velocity is >95% He so a maximum of 5% of the gas 
could condense into dust. Assuming thi s and a different 
gas-to-dust ratio of 100 (Whittet 2003), the amount of 

5 x 10~ b M, 


dust in such a shell would be ~5 x 10~ s M 0 . A different 
gas-to-dust ratio was adopted, since in this scenario the 
envelope is not PN-like. The merger event itself only 
takes a few minute s after which ther e is a rapid expansion 
to supergiant size ( Staff et al.|2012 ). Zhang et al. (20141 


estimate that the new star will expand to ~20U 
500 yr. The estimated velocities of the gas escaping the 
grid in the simulations listed in Table 1 is 600-900 km 
s -1 . This gas would take ~10 4 yr to expand to the size 
of the R CrB shell. These velocities are the same order 


TABLE 2 
Shell Types 


Type 

V exp 

(km s _1 ) 

Time 

(yr) 

Dust Mass 
(M 0 ) 

Total Mass 
(Mq) 

Planetary Nebula 

30 

10 5 

<10“ 3 

<0.3 

WD merger 

600-900 

10 4 

5 x 10~ 6 

10" 2 

RCB star 

400 

10 4 

>10“ 3 



of magnitude as th e winds measured in RCB stars in the 
He I A10830 line (Clayton et al. 20131. This analysis 
assumes that any PIN /common envelope phases of the 
two stars in the binary occurred long before the merger 
but this time period is not known. 


3.4. RCB Phase Mass Loss? 

It is thought that dust forms in pu ffs near the atmo¬ 
sphere of an RCB star (Clayton 1996). We adopted the 
assumption that during a single dust formation event a 
puff forms at 2 R* (R*=85 R 0 ) and subtends a fractional 
solid angle of 0.05, which results in the photosphere of the 
star being obscured a nd a dust mass of ^10~ 8 M 0 (Clay¬ 


ton et al. 


1992 2011). Declines occur when a puff forms 


along the line ot sight, but other puffs are likely form¬ 
ing around the star that do not cause declines. These 
puffs can be detected in the IR. Recent studies pertain¬ 
ing to IR variability of RCB stars has found that the 
covering factor can vary from RCB star to RCB star 
and find an average covering factor of 0.28T0.0 4 for R 


Lambert 

120151). Other studies find a 

tor (Hec 

ht et al. 

1984 

Clayton et 


her co vering tac- 


In addi¬ 
tion,^ acorrelatidn DetweeiT^uIsatioir^hase and the tim¬ 
ing of dust formatio n has been found in several RCB 
stars (Crause et al. 2007) and they typically show reg¬ 
ular o r semi-regular pulsa tion periods in the 40-100 d 
range (Lawson et al. |l990 ). R CrB, itself, does not have 
one regular period but h as shown periods of 40 and 51 d 
(Fernie & Lawson 1993). Therefore, if a dust puff forms 
somewhere around R CrB every 50 days then ~10 -7 M 0 
of dust will form per year around the star. 

There is strong evidence from observations at He I 
A10830 that the dust, once formed, is accelerated quickly 
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by radiation pressure from the star to ^400 km s 


-l 


((Clayton et al. 1992 12003, |2013). The He gas is likely 
;he dust. Thf 


dragged outward by the dust. Therefore, dust produced 
during the RGB phase could fill the R CrB shell in only 
~10 4 yr. In that time, R CrB will produce at least 10 -3 
Mq of dust depending how many puffs form around the 
star. Considering the assumptions made here and those 
made in the radiative transfer modeling, this estim ated 
dust mass is close to the 10~ 2 M 0 estimated by Clay¬ 


ton et al. (2011). The mass of gas in an RCB shell is 
unknown. The gas will be primarily He with little or no 

H. 

The observed far-IR shell around R CrB appears nearly 
spherical. The suggestion that the present day mass loss 
from R CrB is bipolar wou ld not support this morphol¬ 


ogy (Rao & Lambert 


1993). Further polarimetric or in¬ 


terferometric observations are needed to determine the 


morpholo 

gy of the RCB-star dust 

ton et al. 

1997 

Bright et al. 2011 


4. SUMMARY 

The estimated masses of the circumstellar shells in the 
FF/PN, WD merger (DD), and RCB scenarios are sum¬ 
marized in Table 2. If the shell is an old fossil PN in the 
FF scenario, then the gas should be H-rich. In the other 
two cases, it would be H-poor and dominated by He gas. 
The results of the 21-cm observations find no detectable 
H in the R CrB shell. The observations place upper limit 


observed for some PNe are in this range (e.g., Boffi & 

Stanghellini 

1994 

van Hoof & van de Steene 

1999 prte 

radius of R 

JrB shell is very large lor PN shell. The PN 


shells of the FF stars, V605 Aql and Sakurai’s Object, 
are 0.35 and 0.7 pc, respectively. A 4 pc radius PN shell 
would take ~10 5 yr to fill at 30 km s _1 . The gas mass 
loss in a WD merger is ~10 -2 M 0 . Since this gas is 
mostly He, little dust could form in such a shell. 


The gas mass of the R CrB shell is not known, but the 
dust mass has been estimated to be 10~ 2 M 0 (Clayton 


et al. 2011). If it formed during the RCB phase then 
the shell would be filled with He gas, the mass of which 
cannot be measured so the gas-to-dust ratio is unknown. 
Based on the analyses above, we would expect only 10 -3 

10~ 6 M 0 of dust in the PN and WD merger mass loss 
scenarios. Thus the suggestion that the R CrB shell has 
formed from the dust forming during its present RCB 
phase seems most likely since it can form >10 -3 M 0 of 
dust. If true, this model implies that R CrB has been 
an RCB star for ~10 4 yr to have produced a 4 pc radius 
shell. More sensitive 21-cnr observations of RCB star 
shells are needed to place more stringent constraints on 
their H masses. 
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